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Abstract. We consider the fluctuations of the free energy of positive temperature directed poly- 
mers in thin rectangles (N,N a ), a < 3/14. For general weight distributions with finite fourth 
moment we prove that the distribution of these fluctuations converges as N goes to infinity to the 
GUE Tracy- Widom distribution. 

1. Introduction 

We prove a statement of KPZ universality for directed polymers in random media in thin rect- 
angles (N,N a ) for a < 3/14. We assume general weight distributions with at least a finite fourth 
moment. Our main result is the positive temperature analogue of Corollary 1.1 of [I], Theorem 1 
of [3], and the result of [18] . where a similar universality for last passage percolation was proved. 

The model we consider is defined as follows. 

Definition 1.1. Let Wij, i,j G N, be a family of i.i.d. random variables with E[Wii] = 0, 
E[(Wn) 2 ] = 1 and E[(W n ) 4 ] < oo. For each j,k G N we define S j {k) = Ya=i w H with the 
convention 5 J (0) = 0. The partition function for the discrete directed random polymer from (1, 1) 
to (JV, n) at inverse temperature (3 > is defined as 

ZnM= E exp^f;(^^)-^fe_i-l))Y (1) 

l=i <i 1 <...<i n =N ^ j=l ' 

Remark 1.2. The hypotheses that E[Wu] = 0, E[(lfn) 2 ] = 1 are not restrictive. Indeed, a non- 
zero mean only changes Z^^ n {(5) by a (deterministic) multiplicative constant while a different second 
moment is just a rescaling of (3. Also note that one can rewrite Z^ tn {f3) as the sum over all up/right 
lattice paths ir from (1, 1) to (JV,n) of the Boltzmann weights exp(/3 j)67r ^y)- 

The main result of the paper is the following: 

Theorem 1.3. For all a G (0, ^) and P>0,ifn= [N a \ then 

JV->oo V /3JV3-6 / 

for each fixed r£l, where Fque is the GUE Tracy-Widom distribution function [19J. 



Remark 1.4. If we assume that the Wij have all moments finite, the above result holds for all 
a < |. See Remark 11.71 below. 

The thinness assumption, n = [N a \ with a G (0, A), is only a technical condition. It is antic- 
ipated that a universality result should also hold for thicker rectangles (for all < a < 1), after 



i 
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an appropriate modification of the centering and scaling of log Zj\r n (/3)- The modification should 
depend on the distribution of weights; however, there is presently no prediction for the exact de- 
pendence for the full rectangle case when n = AU, For a < 1, at least the leading order term of the 
centering was already known. Based on the earlier works of |9] [16] and estimates similar to some 
of those in the proof of Proposition 11,61 Moreno Flores \12\ Theorem (8.8)] showed a law of large 
numbers for the free energy by proving that with probability one: 



We prove Theorem 11.31 by employing the Skorohod embedding theorem to couple the discrete 
directed polymers to the semi-discrete directed polymer of 0'Connell-Yor[15j. The asymptotics of 
the free energy of the semi-discrete polymer is then evaluated via a steepest descent analysis of the 
Fredholm determinant formula of [H Theorem 5.2.10]. 

This coupling approach is an adaptation of the previous works on directed last passage percolation 
by [I] and [3], which corresponds to the case when /? = oo. We note that the same thinness 
condition a < 3/14 was assumed in these works under the finite fourth moment condition. Let 
Ln,tl = hm^^oo /3 _1 log Z]\r jn ((3) denote the last passage percolation time. Then it was shown that 



This is consistent with our main result if we formally take the limit /3 — > oo and replace /3 1 log Z/v, n (/3) 
by L N>n , 

When we assume that the first p moments are finite for p > 2, then it was shown in [3] that (J3j) 
holds fora<|(| — |). This is obtained by using the Komlos-Major-Tusnady theorem instead of 
the Skorohod embedding theorem in the couple argument. Note that the upper bound is | when 
p = oo, which is in agreement with the result in this paper stated in Remark 1 1 . 4 1 ab ove . We remark 
that there is a very different proof of (|3|) by [18] assuming that p = 3. This proof is based on a 
concentration inequality of [6]. The condition on a is same asa<|(^ — |) with p = 3. It is an 
interesting question to find this approach can also be adapted to the polymer case. 

The semi-discrete directed random polymer of O'Connell-Yor [15] is defined as follows. 

Definition 1.5. Fix t > and n £ N. To each ordered set = to < t\ < ■ • • < i n _i < t n = t, we 
associate an up/right semi-discrete path <j) = {(x,i) : < x < t{\ in [0, t] x {l,...,n}. Given 
a family of independent standard one dimensional Brownian motions B*, 1 < i < n, we define the 
energy of an up/right path (ft as 



There is presently only one known exactly solvable positive temperature discrete polymer. It has weights which are 
distributed as the logarithm of inverse Gamma random variables. Introduced by Seppalainen [IT], this polymer has 
an explicit product invariant measure which allowed him to compute the law of large numbers for the free energy as 
well as a tight upper bound on the free energy fluctuation scaling exponent for this model. Due to the tropical RSK 
correspondence [Jj, this polymer also fits into the hierarchy of solvable models studied in [4]. Though it has not yet 
been done, this should enable a rigorous proof of the GUE Tracy- Widom limit theorem for any size rectangle as long 
as both n and iV go to infinity. Due to the work of [14], another model which fits into the hierarchy of [4] is the 
semi-discrete directed polymer model, introduced by O'Connell-Yor [15] (see discussion below Proposition 11.80 . 




-^logZ 7Vin (/3) = 2/3. 

2 + 2 




(3) 



E(<j>) = B\t!) + (B 2 (t 2 ) - B 2 (ti)) + ■■■ + (B n (t n ) - B n (i n -i)). 
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The partition function for the semi-discrete directed random polymer at inverse temperature (3 > 
is then defined as 



?' 



{fi) = J exp (pEtf)) d<f> 



0=to<tl<-<tn-l<tn=t 

where d<\> denotes the Lebesgue measure in over the variables t\, . . . , t n -\. 



Theorem 11.31 is a consequence of the following two results which we prove in Section [2j The first 
is a coupling result of Z^ jn {(3) and Z^{(3). Here we do not assume anything on the relationship 
between n and N. 

Proposition 1.6. Given 5 > 0, there exist positive constants C = C(5,(3) and c = c(/3) such that 



for all n,N G N and a > 0. 

Remark 1.7. If we assume that the Wy have all moments finite, then we can couple the discrete 
and semi-discrete polymers using the Komlos-Major-Tusnady theorem and replace the term N± +s 
to N s . For further detail, see Remark 12.51 It is easy to check that this improvement is enough to 
prove Theorem 11.31 for a < |. 

The second result is on the asymptotics for Z^((3). 

Proposition 1.8. For each fixed a G (0, 1) and /3 > 0, 

/logZ^ J (l)-2/% 1 - K \ 

for every fixed ret, where k = and [i = ^p- 

Note that for a £ (0, 1) we have that k G (0, 1/2) and /i G (1/3, 1/2). For a = 1 an analog to 
the above result with fi = 1/3 is proved in ^ Theorem 5.2.12] and Theorem 2.1]. In this case 
the form of the law of the large numbers term which must be subtracted is more involved. Our 
proof of Proposition 11.81 is an adaptation of the analysis of [H [5]. On the other hand, for a = 0, if 
the number of rows is kept fixed at n then the semi-discrete polymer free energy converges (after 
fjL = 1/2 scaling) to a semi-discrete last passage percolation time. It is known that this last passage 
time is exactly distributed as the largest eigenvalue of an n x n GUE random matrix (see [21 [10] ; 
see also Theorem 1.1 of |14j). Thus, when a = 0, if we take t to infinity first and then let n tend to 
infinity, we also obtain the -Fgue distribution as the limiting fluctuation. 

We can now prove Theorem 11.31 

Proof of Theorem \1.3[ We apply Proposition II .81 with t = N and note that Brownian scaling implies 
that Z™(j3) is equal in law to /3 -2 ( n_1 ) Z% t (X) ■ This reduces the proof of Theorem 11.31 to the claim 
that 



logZ N>n ((3)-logZ%((3)-2(N a - 1) log /3 



converges to zero in probability. 
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We may neglect the term 2(N a — 1) log /3 due to the fact that, 

N a < JVs-f for all a < 3/7 

as iV~ goes to infinity. 

1 a 

For any fixed e > 0, we apply Proposition 11.61 with a = eN^^e 
5 > there exist positive constants c and C such that 



and n = [N c 



(4) 

to see that for all 



logZjv.nOS)- logZftOS) 



> e 



where ?y(iV) = cN~ 2a (a - log(iV a !)) 2 . By Stirling's approximation, log(iV a !) » iV a logiV a - TV". 
Using (JH) again, we find that a ^> log(iV a !) for N large, and hence rj(N) > c'iV 1 ^ for some other 
constant d depending on c and e. 

As iV goes to infinity, the second term converges to zero if a < | and the first term converges 

_3_ 

14 

Theorem 11.31 



to zero if a < A by choosing 5 to be small enough (but fixed in N). This completes the proof of 

□ 
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2. Proofs 

We prove Proposition 11.61 and Proposition 11.81 



2.1. Proof of Proposition 11.61 We compare the discrete polymer of general weight distribution 
(satisfying the conditions of Definition II. ip with the semi-discrete polymer by coupling them using 
the Skorohod embedding theorem. 

For much of the proof it will be convenient to rescale our processes so as to be functions in 
C([0, l],M n ), the space of continuous functions from [0,1] to M ra . We will denote such rescaled 
processes with an overbar to distinguish them from the unsealed versions. For two such functions 
/ = (/!,...,/«) and g = {g\...,g n ) in C([0,l],M n ) define the metric 

n 

d{f\g) = Y J \\f j -g'\U (5) 

i=i 

where 

11/^-^1100= sup \p(t)-gi(t)\. (6) 
te[o,i] 

Set 



F N (f) = log 



rt 



£ exp^C/^)-/^-^)) 

■ff=to<tl<—<t n -l<t„=l 



(7) 



for N £ N. A basic property of F/v is the following: 



Lemma 2.1. The function Fn ■ C([0, l],M n ) —> K defined in ((7]) is Lipchitz continuous with Lipchitz 
constant less than or equal to 2/3 with respect to the metric ©. 
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Proof. Set A(f) = e FN ^\ Triangle inequality implies that A(f) < A{g)e 2 ^ d ^ ,g "> . Therefore, assum- 



ing without loss of generality that > 1) we have \F(f) — F(g)\ = log 4^4 < 2(3d(f,g). □ 



Ma) 



Mf) 



Ms) 



The next lemma gives us a coupling between the random walks S 3 (k) in Definition 11.11 and the 
Brownian motions B 3 {t). This is the Skorohod embedding theorem (whose proof can essentially be 
found in [81 Theorem 8.6.1]) applied to each walk S 3 (-). 



with respect to the Brownian filtration of B 3 ) such that E[r^] = 1, E[(r^) 2 ] < oo, and 



Lemma 2.2. For each j, there exist i.i.d. stopping times tI, k = 1,2,..., (which are measurable 

i\ = i> n(ri) 

Bi( T { + ---+T 3 k ) = S\k), fc = 1, 2 

Define 

= B 3 {r{ + ■■■ + r{ m ) + (Nt - [Nt\)[B 3 (r( + ■■■ + rf m] ) - B 3 {r{ + • • • + r( m )) (8) 

to be the piecewise linear interpolation of the Brownian motion sampled at times t\ + • • • + rl for 
I < k < N , rescaled so as to be a function in C([0, 1],R). 

Consider the continuous vector valued function SV(i) ■= (Sjj(t), . . . , S^(t)) G C([0, l],M n ). From 
LemmaEa S j N (±) = S 3 N (k) for 1 < k < N. Therefore, from CD), and ©, one finds that 

logZ N>n (p) = F N (S N ). 

We also define the rescaling of the Brownian motions by 

B N (t) = (B}r(Nt),...,B%(Nt)). 

The proof of Proposition 11.61 is now obtained by estimating probabilitistic bounds on F/v (Sn) — 
Fn(Bn) and Fn(Bn) — logZ]^(/3). The relevant estimates are obtained in Lemma [2.31 and in 
Lemma 12.61 respectively. 

Lemma 2.3. For any given 5 > 0, there exist a positive constant C such that 

f(\f n (s n )-f n (b n )\ > a \ < ^m+ 5 

for all n,N G N and a > 0. 

Proof. From Lemma 12.11 the Markov inequality, and the union bound 

\F N (S N ) - F N (B N )\ >a) <v(d(S N ,B N )>^\ < ^E||5^ - B^loo- (9) 



Define the event 



An(u, A) := <j sup 

Ki<N 



> uN 1 - x 



for u > and A G (0, 1). From Doob's inequality ([8l Theorem 5.4.2 and Example 5.4.1]) and the 
definition of we have 

P (^,A))<^i^. (10) 
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We now fix p > and write 



S l N -B 1 N \\ 00 >uN p \ du 



<N P + N P 



\Sat — B 



N ~ -"AT I loo 



> uN p du. 



(11) 



In order to estimate the last integral, we fix A € (0, 1) and divide the sample space into A^{u, A) 
and An(u,X) c . From (flOl) we find that 



N p 

It remains to estimate 

/•OC 

N p 



S l N - S^Hoo > uN p , A N {u, A) ) du< Var (t{)N p 



-1+2A 



(12) 



> uN p , A N (u,X) c du. 



\S]\f — BjyWoo 

For this purpose, we note the following result of P. Levy on the modulus of continuity of a Brownian 
motion. (See the proof of [81 Theorem 8.4.2].) 



Lemma 2.4. There are positive constants K\,Ki such that 



sup \B 1 (s)-B 1 (t)\ > x < — e 



s,te[o,i] 

\t-s\<r 



(13) 



for all r E (0, 1) and x > 0. 

By ([8]) and the triangle inequality, 

\S 1 N (t)-B 1 N (t)\\ 00 >uN p , A(u,XY 

ip S 1 (T 1 1 + ---+r; iVtJ )- J B 1 (iVt) >\uN p , A(u,\) c ) 

3,1] 2 / 

^Vi + • • • + r r Vl ) " # V! + • • • + r L Vj )| > \uN p , A(u, A) c ) . 



sup 

te[o, 



(14) 



+ P sup 
\te[o,i] 

On A(X, u) c we have l^ 1 + • • • + T^ m] - [Nt\ \ < uN x - x for all t E [0, 1]. Since A E (0, 1), we have 
uN 1 - x > iV 1 " A > 1 for all u > 1 and N > 1. Hence on A) c we have for all i E [0, 1] and N > 1 

l^ 1 + . . . + r L V tJ - ATt] < 2uiV 1 - A . 



(15) 



Thus, the first term in the right side of (|14p is bounded above by 



( sup ^(t) - B 1 (s)\ > \uN p \ 



,t£[0,N+uN 
\t~s\<2uN 1 ~ x 

Similarly, the second term in the right side of (|14p is bounded above by 



sup \B 1 (t)-B 1 (s)\ > ^unA. 



(16) 



s,te[0,N+uN 
\t-s\<3uN 1 - x 
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We must consider s, t € [0, N+uN 1 A ] because on the event A(u, A) c , (Tj L +. . .+Tnv£i ) < N+uN 1 x . 

We wish to estimate the above probabilities using Lemma 12.41 To do that we must rescale time 
by the factor N + uN 1 ~ x . Thus define t = (N + uN^yH. Then B 1 ^) = {N + uN 1 ^) 1 / 2 ^ 1 ^) 
where B 1 is another standard Brownian motion. This rescaling reduces (|15p to the left-hand side 
of with 

2uiV 1 - A _ U JV 

r ~ iV + uiV 1 ^' X ~ 2 V ^VT^iV T ^' 
Applying Lemma 12.41 shows that 



Likewise one sees that 



3 V u 



Integrating the sum of these bounds over the interval u G (1, oo) and multiplying N p , we find 
that for some other positive constants K% and K4, 



NP J™F^\\Slf -BjfWoo > uN p , A N (u,X) c ^j du < e~ K ^ P ~ 



l+A 



for all JV > 1. 

Combining with (jlip and (|12p . we conclude that 

E||^ - B^lloo < NP + Var (rf)^-^^ + K^N 1 -^*^^ . (17) 

Recall that Var (rj 1 ) < 00 by Lemma [2721 Now choosing A = ^ and p = \ + S with (5 > 0, Lemma 
follows from @ and flU]). □ 



Remark 2.5. If we assume finite exponential moments for the weight distributions then we can 
better approximate the simple random walk by a Brownian motion using the dyadic approximation 
of Komlos-Major-Tusnady. Indeed, using (111 Theorem 7.1.1], one sees that in this case for any 
6 > there exists a constant C > such that 

Ep^-B^Hoo < CN 5 . 

Then from ([9]) the term N± +s in Lemma 12.31 is improved to N s . This in turn implies the same 
change in the first term of the bound in Proposition 11.61 It is easy to check that with this improved 
bound the proof of Theorem 11.31 goes through under the weaker assumption that a < | . 

We now compare the partition function of the semi-discrete polymer with Fn(Bn)- 

Lemma 2.6. There are positive constants K\ and K% such that 

/ _ \ A' 2 (a-log(n!)) 2 

F l\F N (B N ) - log Z^(/3)\> a) < KmNe ^ . 
for all h,JVeN and a > 0. 
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Proof. The proof follows the same strategy of Lemma 12.31 We will use Lemma 12.41 For u > 0, we 
define the event 

S(u) = ( sup \B j (t) - B j {s)\ > u 

0<|t'-s'|<l 

By Lemma 12.41 Brownian scaling, and the union bound, 

P(S») < K in Ne- K2U \ (18) 



Let us now observe that we can write 

ZW)= £ A(t ,...,t n ) 



l=t <t 1 <-<t n - 1 <t n =N 

Uez 



(19) 



where 

A(t ,...,t n ):= f e^d<t>. 

J D(t ,...,t„) 

In order to define the domain D(to, . . . , t n ) over which (j) is integrated, recall that (f> can be specified 
by an ordered set = so < si < • • • < s n -i < s n = N. Then D(to, . . . ,t n ) is the set of <fi whose 
associated ordered set also satisfies Sj G [U — 1, ti) for 1 < i < n — 1. 
On the event S(u) c , we have 

Integrating the above inequality over <fi G Z)(to, • • • , t n ) and noting that the Lebesgue measure of 
D(tQ, . . . ,t n ) is bounded in [(rt!) , 1], we find that on S(u) c , 

r rl n-l e /9E"=i(B :, '(%)--B 3 'fe-i-l))-2 J 8n« < . . . )in ) < e /3E;=i(fi J fe)-B J (t J -i-l))+2/3n« > 

Hence on S(u) c , 

F N (B N ) - 2f3nu - log(n!) < logZ^(/3) < Fjv(5jv) + 2(3nu. 
Thus from (1181) we find that 



P( \F N (B N )- log Z%(p)\ > 2/3nn + log(ra!) ) < P(5(«)) < i^iniVe-^ 2 " 



for all u > and n, N G N. Setting a = 2/3nn + log(n!), we obtain the claimed result of this 
lemma. □ 

The proof of Proposition 11.61 follows immediate by combining the above Lemmas (I2.3P and (12.6p . 

2.2. Proof of Proposition 11.81 This proof is based off of the proof of Theorem 4.1.46 of [4] 
which establishes the GUE limiting distribution for the semi-discrete polymer with a constant ratio 
between n and t. Both that result and the present result rely upon Theorem 5.2.10 in [4 J which 
(specializing that result to the case of a± = ■ • • = ajv = 0) yields: 

Theorem 2.7. Fix n > 1 and r > 0, and < S < 1. Then, 



E 



e -«S?(i) 



det(/ + if u ) 
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where K u : L 2 {Cq) — > L 2 {Cq) is the operator defined by the kernel 

Ku ( v v >) = J_ ds r(- s )T(l + s) - [ ' r — 

and Co is a positively oriented contour containing and such that \v — v'\ < 5 for all v,v' € Cq. 

The above result for the Laplace transform leads to the desired asymptotic probability distribution 
by virtue of the following result which is given as Lemma 4.1.38 of [3J. 

Lemma 2.8. Consider a one-parameter family of functions {ft}t>o mapping M. — > [0,1] such that 
for each t, ft(x) is strictly decreasing in x with a limit ofl at x = —oo and at x = oo, and for each 
5 > 0, on M\ [—5, 5] ft converges uniformly to l(x < 0). Define the r -shift of ft as fl{x) = ft(x — r). 
Suppose that a one-parameter family of of random variables X t satisfies, for each r£l, 

lim E[f[(X t )]=p(r) 

t— ¥00 

for a continuous cumulative distribution function p(r) . Then Xt converges weakly in distribution to 
a random variable X which is distributed according to F(X < r) =p(r). 

Let fi > and n > be defined as in the hypothesis of Proposition 11.81 Consider the functions 
ft{x) = e~ e . Observe that this family of functions meets the criteria of Lemma 12.81 By Lemma 
I2.8[ if for each r £ R we can prove that 



lim E 

t— too 



rr 
■It 



log 4^(1) -2/3* 1 -" 



then it will follow that 

/logZj*" J (l) - 2pt 1 - K \ 
lim P( P t \' < r ) = F GUE (/rV). 

Observe that if we define 

u = e- 2 ^- K - rt \ (20) 

then 



f 

Hence, in view of Theorem 12.71 our proof reduces now to proving that for u as in (|20p . n = \t a \ 
and t = f3 2 t 

lim det(J + K u ) = F GUE (/rV) (21) 

t— >oo 

for each fixed r£l, 

We prove ([21]) using a steepest descent analysis. This analysis follows exactly the approach 
employed in the proof of Theorem 4.1.46 in [3]. As such we only include the critical point analysis 
here presently. The necessary contour manipulations and tail bounds can readily be found in [3J 
and easily adapted. In the below we replace n = \t a \ by n = t a to make the notations simpler. It 
is straightforward to check that the error coming from this change is negligible 

Let us consider the kernel K u (v,v'). We scale the kernel by setting v = t~ K v and v' = t~ K v' and 
define K u (v,v') = K u (v,v')t~ K . Then det(I + K u ) = det(I + K u ). By changing the variables as 
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s = t K (C — v) and v = t K v in the formula of the kernel, we have 

K u (v,v') = ^- [ , exp {t<* (G{v) - G(C)) + rtf (v - C) } J^j (22) 

where 

G(-z) = lo g r(^ K z) - l3 2 z 2 /2 + 2/3 z 

and we have also replaced T(— s)r(l + s) = tt/ sin(— its). Here the contour for £ in the integral may 
be taken as v + 5 + lR for any fixed 5 > and the operator K u may be defined on L 2 (C*o) where 
Co is a positively oriented contour containing and such that \v — v'\ < 5 for all v,v' £ Co. The 
problem is now prime for steepest descent analysis of the integral defining the kernel above. 

The idea of steepest descent is to find critical points for the argument of the function in the 
exponential, and then to deform contours so as to go close to the critical point. The contours 
should be engineered such that away from the critical point, the real part of the function in the 
exponential decays and hence as t gets large, has negligible contribution. This then justifies localizing 
and rescaling the integration around the critical point. The order of the first non-zero derivative 
(here third order) determines the rescaling in t which in turn corresponds with the scale of the 
fluctuations in the problem we are solving. It is exactly this third order nature that accounts for 
the emergence of Airy functions and hence the Tracy Widom (GUE) distribution. 

The Digamma function is defined as ^f(z) = [logr]'(z) and, along with its first two derivatives, 
have the following expansion for \z\ small: 

*(*) = — + O(l), V(z) = \ + 0(1), <f"(z) = -4 + 0(1). 

z Z A z 6 

Let us then record the first three derivatives of G along with their large t expansions (afforded by 
the above expansions) 

G'(v)= ^(t- K i)t- K - (3 2 v + 2p = — - p 2 v + 2j3 + 0(t- K ) 

v 

G"{v)= ^'{t- K v)t- 2K - /3 2 = 4t - f + 0(t- 2K ) 

v z 

G"'(v) = ^"(t- K v)t- 3K = -4 + 0(t- 3K ). 

The critical point is the value of v = v c at which G'(v c ) = 0. However, it suffices to choose v c 
such that 

(3 2 v c + 2/3 = 

since then G'(v c ) = 0(t~ K ). Solving the above gives v c = This choice of v c implies that 

G'(v c ) = 0(t- K ) 

G"(v c ) = 0(t- 2K ) 

G"'{v c ) = -2/3 3 + 0{t- 3K ). 

Note that by taking 5 small enough, it is possible to deform the contour of ([22]) to pass the critical 
point C = P ■ 
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Therefore we may make the final change of variables to expand around v c by setting v = v c +t~^v. 
From Taylor approximation and the above bounds on the derivatives of G we find that 

G(v) = G(v c )-^-t~ a (vf + lot 

where lot denotes lower order terms in t. We also set £ = v c + t~ 3£. We may conclude from the 
above argument that (note that the constants G(v c ) cancel) 

t a (G(v) - G(0) = -^-(v) 3 + ^(C) 3 + o(l). 



3 v ' 3 

Now we can turn to the remaining parts of the kernel. Note that 

rt%(v — C) = r {v — C)- 
As for the other two terms in K u we have 

-Kt' K 1 d( d( 



r 



sin(7rt K (v — Q) v — ( ( — v' Q — v' 

where the t~~s comes from the Jacobian associated with the change of variables from v to v. 

Putting together the pieces we have that (modulo the necessary tail estimates and uniformity 
which can be obtained as in the proof of Theorem 4.1.46 in [4]) 

lim det(I + K u ) = det(I + K r ) 

t—^oo 

where 

K r : L 2 (C V ) -> L 2 (C V ) 

for C v a contour given by rays from the origin at angles ±27r/3 oriented to have increasing imaginary 
part. The operator K r is defined in terms of its integral kernel 

1 f e -^« 3 +™ d( 



c< e-ff+'i (8 - OK - •>') 
where Cq is a contour given by rays from d (for some d > 0) at angles ±7r/3 oriented to have 
increasing imaginary part. But it is known that (see Proof of Theorem 4.1.46 in [3]) 

det(I + K r ) = F GUE (/5-V). 

This proves (|21|) and hence completes the proof of the proposition. 
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